Universality of Level Spacing Distributions in Classical Chaos 
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We suggest that random matrix theory apphed to a classical action matrix can be used in classical physics 
to distinguish chaotic from non-chaotic behavior. We consider the 2-D stadium billiard system as well as the 
2-D anharmonic and harmonic oscillator. By unfolding of the spectrum of such matrix we compute the level 
spacing distribution, the spectral auto-correlation and spectral rigidity. We observe Poissonian behavior in the 
integrable case and Wignerian behavior in the chaotic case. We present numerical evidence that the action matrix 
of the stadium billiard displays GOE behavior and give an explanation for it. The findings present evidence for 
universality of level fluctuations - known from quantum chaos - also to hold in classical physics. 
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Introduction. Experiments and computer simulations 
have shown quantum systems with a fully chaotic classi- 
cal counter part possess universality properties IH 0, S 
Ulstl- For example, the spectra originating from differ- 
ent heavy nuclei give the same Wignerian energy level 
spacing distribution [6]. About 1950 Mehta and Wigner 
showed that random matrices of Gaussian orthogonal 
ensembles (GOE) generate a Wigner- type level spac- 
ing distribution |1J. The celebrated Bohigas-Giannoni- 
Schmit (BGS) conjecture |7|J| states that in time-reversal 
invariant quantum systems with fully chaotic classical 
counterpart, the energy level spacing distribution is the 
same as that obtained from random matrices of a certain 
symmetry (Gaussian orthogonal ensembles GOE), re- 
sulting in a Wignerian distribution. This paper is about 
classical chaos occuring widely in nature, for example 
in astro physics, meteorology and dynamics of the at- 
mosphere, fluid and ocean dynamics, climate change, 
chemical reactions, biology, physiology, neuroscience, 
or medicine. Traditionally, classical chaos is described 
by tools of nonlinear dynamics like Lyapunov expo- 
nents, Kolmogorov-Sinai entropy and phase space por- 
traits (Poincare sections). In this work we present evi- 
dence that fully chaotic classical systems show univer- 



sality behavior, visible in level spacing fluctuations ob- 
tained from a classical action matrix. We show that ran- 
dom matrix theory can be used in classical systems to 
distinguish chaotic from integrable systems. In the case 
of the 2-D billiard system, Lorentz gas and anharmonic 
oscillator we find GOE-type behavior of action matrix 
elements and Wignerian behavior in the level spacing 
distribution. Moreover, we provide an explanation why 
those action matrix elements display GOE-type behav- 
ior This holds when the system is both, classically 
chaotic and diffusive (random), like the Lorentz gas and 
the reason is due to the Central Limit Theorem. We dis- 
cuss implications on ergodicity. 
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FIG. 1: Scheme of 2-D stadium billiard system with trajecto- 
ries going from x, to xj located on the wall of a quarter billiard. 
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Classical systems: action matrix and level spacing 
distribution. In quantum systems a discrete level spac- 
ing distribution occurs, if the potential supports a bound 
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FIG. 3: Histrogram of action matrix elements. Top: stadium 
billiard corresponding to 12 collisions. Bottom: diffusion 
(random walk) corresponding to 1000 collisions. A Gaussian 
behavior is predicted by the Central Limit Theorem. 



FIG. 2: Stadium billiard system in 2-D. Level spacing distri- 
bution (top) from action matrix £ and its spectral rigidity A3 
(bottom). 

state spectrum, like in the hydrogen atom. The classical 
Hamiltonian function H{x,p,t) yields a continuous en- 
ergy spectrum, and hence is not suitable to play a role 
like the Hamiltonian of the quantum system. However, 
the Hamiltonian function H{x,p,t) is closely related to 
the Lagrange function L{x,x,t). In order to construct 
a finite matrix, we consider the classical action E evalu- 
ated along classical trajectories between discrete bound- 
ary points, 

E = S[x"-"j] = L{x{t),x{t),t) 1,=,™; . (1) 

Jo 

We choose a finite set of discrete points {xi,x2, . . . ,xn}. 
For any pair of boundary points jc/,xy we compute the 
corresponding classical trajectory Xx"Jj (we assume for 
the moment that to each pair of boundary points corre- 
sponds a unique classical trajectory, the case of mul- 
tiple trajectories will be discussed below). Then the 
action matrix is built from matrix elements given by 



E,y = S[xx"xj]- We then compute the spectrum of eigen- 
values a, of the action matrix E, ; and carry out a statisti- 
cal level spacing analysis yl |4], |5|] . Here we want to ex- 
plore the hypothesis that such action matrix for chaotic 
classical systems plays a role analoguous to the Hamil- 
tonian matrix for chaotic quantum systems. 

Billiard System. Let us consider as example the sta- 
dium billiard in 2-D. The billiard is symmetric under 
mirror operation about x- and y-axes at the center. In or- 
der not to mix different symmetry classes, we consider 
the billiard system (see Fig.ITl) with boundary points 
located only in one quarter of the billiard. Trajectories 
start from and arrive at points Xi and xj, respectively, lo- 
cated at the billiard wall. The rule of dynamics is free 
motion in the interior region and perfectly elastic specu- 
lar collision at the wall. Such a classical billiard system 
is known to be fully chaotic For a given pair of 
boundary points, there is an infinite number of trajec- 
tories, connecting these points. We have discarded tra- 
jectories of small angle collisions (staying close to the 
wall). We have classified the trajectories according to 
number of rebounds Nreb from the stadium walls. One 
expects chaoticity to increase with Nreb- However, from 
the computational point of view, Nreb can not be too 



3 



large, because each rebound from the curved part of the 
wall decreases the numerical precision roughly by one 
order of magnitude. The results presented below corre- 
spond to Nreh of the Order of 4 to 6, which turned out 
to be sufficient to represent the chaotic behavior We 
have taken care that numerical values of action matrix 
elements have a relative error below 10^^. For the com- 
putation of trajectories we used a symplectic algorithm. 
The numerical results for the level spacing distribution 
P{s) of action eigenvalues and the spectral rigidity A3 
are shown in Fig.Q. They represent results of super- 
imposition of spectra corresponding to different shapes 
of the billiard. We considered shapes parametrized by 
horizontal over vertical diameter of the stadium (vary- 
ing from 2.2 by increments of 0.1 to 3.2). For each bil- 
liard shape we considered 50 boundary points evenly 
distributed on the wall of the quarter billiard and com- 
puted trajectories for all pairs of boundary points. To 
unfold the spectra we used the technique of Gaussian 
broadening isl The results are consistent with a 
Wigner distribution. 

GOE behavior of action matrix elements in billiard 
system. In quantum chaos the BGS conjecture postu- 
lates equivalence, with respect to energy level fluctua- 
tions, of matrix elements from the (time-reversal invari- 
ant) quantum mechanical Hamiltonian and GOE ran- 
dom matrix elements. Here we show that classical ac- 
tion matrix elements of the stadium billiard behave like 
GOE matrices. We have investigated the statistical be- 
havior of action matrix elements for the original stadium 
billiard. The action matrix elements were computed for 
long trajectories (12 rebounds). In GOE all off-diagonal 
matrix elements (and likewise all diagonal matrix ele- 
ments) obey the same distribution. Thus, for statistical 
purpose, one can draw values from different matrix el- 
ements. Likewise in case of the stadium billiard, we 
computed the action for an ensemble of different ma- 
trix elements and compiled a histogram. The result is 
shown in Fig. 131, and gives approximately a Gaussian 
(however, with non-zero mean). 

In order to explain such behavior, we present here an ar- 
gument (not a proof) that this is due to the Central Limit 
Theorem. Let us consider the full billiard. The particle 
(representing the billiard ball) alternates between colli- 
sions with the semi-circle walls and free motion (colli- 
sions with straight walls are inessential because the sin- 
gle billiard is equivalent to an infinite number of bil- 
liard copies attached side-by-side at the straight walls 
and then removing the straight walls). Such system is 



very similar to the system of a particle colliding with 
fixed disks: the so-called 2-D Lorentz gas. Here the 
particle alternates between free motion and collisions 
with circular disks. Furthermore, the Lorentz-gas model 
of equal disks located on a rectangular grid is equiva- 
lent to the Sinai billiard. The Lorentz gas model has 
been shown to have chaotic character (positive average 
Lyapunov exponent) and diffusion character (finite dif- 
fusion constant) [lOfl. The Lorentz gas model has a dif- 
ferent character in the case of dense packing of disks, 
which, however, we do not consider here. In the Lorentz 
gas model, in the limit of many collisions between parti- 
cle and disks, the particle can be viewed as an isotropic 
random walker. Assuming that a particle moves in fixed 
time steps A f , at each step a distance A x drawn from a 
uniform probability distribution over the interval [0, 5], 
and travels in a direction of an angle drawn from a uni- 
form probability distribution over the interval [0, 2jt], 
the probability density function for the final position of 
the particle obeys the diffusion equation dtP— DV^p for 
sufficiently large times. Accordingly, the variance in the 
final position of the particle grows linearly in time, i.e. 
(x^) ~D t where D is the diffusion constant. The action 
of the random walker going from xq to x/ is 

^ 2 Af ' ^ ^ 

Thus Eg is a sum of random variables, all drawn from 
the same distribution. The central limit theorem implies 
that is a random variable obeying a Gaussian dis- 
tribution (in the limit 0°). This is confirmed by a 
numerical simulation shown in Fig. [S) . Note that the av- 
erage total action for the isotropic random walker is non 
zero and depends on its mass and mean free path, i.e. 
the average distance travelled in one time step. Since 
the action Eq of a random walker is distributed accord- 
ing to a Gaussian distribution and a particle in a Lorentz 
gas system presents diffusive behavior, we expect that 
the action Eq of a particle in a Lorentz gas system is 
distributed according to a Gaussian distribution. On 
the other hand the Lorentz gas model is a determinis- 
tic classical Hamiltonian system obeying time-reversal 
symmetry. Hence the action matrix elements Eq form 
a GOE-type ensemble (with non-zero mean value and 
variance depending on the parameter A^). Because GOE 
random matrices yield a Wignerian level spacing distri- 
bution yj], we expect the diluted Lorentz gas model, in 
the regime of many colUsions, to yield a Wignerian level 



4 



spacing distribution from the classical action matrix. 

Harmonic and Anharmonic Oscillator. We have also 
considered the 2-D harmonic oscillator (mass ni and fre- 
quency (o) and also the 2-D anharmonic oscillator (with 
coupling 'hp-y'^). We have chosen the following param- 
eters: m = (0 = 1 and A- = 8. By computing the Poincare 
sections (which depend on the energy) we have veri- 
fied that the system is classically almost fully chaotic. 
We superimposed spectra corresponding to variation of 
0) and X, such that co/A, = 1 /8, and X taken from the 
interval [8, 30]. We computed the action matrix E,y, 
unfolded the spectrum and computed the level spacing 
distribution, as well as the spectral rigidity A3. The 
results for the anharmonic oscillator (not shown here) 
are consistent with a Wigner distribution of level spac- 
ings and also consistent with GOE behavior in spec- 
tral rigidity. In contrast, the harmonic oscillator (inte- 
grable system) yields a level spacing distribution and 
spectral rigidity of Poissonian nature (results not shown 
here). In smoothing the raw spectrum and unfolding the 
spectrum, we found it useful to tune the free parameter 
such that to reproduce the auto-correlation coefficient 
(C = for Poissonian distribution and C = —0.27 
for Wigner distribution HI). Meeting this condition 
helped to fine-tune the level spacing distribution to be- 
come closer to Poissonian, or Wignerian, respectively. 

Discussion. We have suggested to extend random ma- 
trix theory, used in chaotic quantum systems, to classi- 
cally chaotic systems. We found for the classical har- 
monic oscillator in 1-D and 2-D (integrable system) a 
Poissonian action level spacing distribution. In contrast, 
for the 2-D classical anharmonic oscillator in the chaotic 
regime as well as for the classical stadium billiard we 
found a Wignerian action level-spacing distribution. 
From these findings we are led to propose the follow- 
ing conjecture: Spectra from action matrices of time- 
reversal invariant classical K-systems show the same 
statistical fluctuation properties as those predicted by 
GOE matrices. If this holds, it means that universal 
laws hold in fully chaotic classical systems. We have 
presented a plausibility argument, based on the central 
limit theorem, why action matrix elements of the sta- 
dium billiard should obey GOE-type laws. 
The diffusion process (random process) and determinis- 
tic chaotic motion of the Lorentz gas are different with 
respect to memory: In the random process memory is 
lost in one collision step, while in the chaotic system 
memory is eradicated exponentially with the number 
of collisions (given a finite resolution). However, the 



chaotic Lorentz-gas system, after sufficiently many col- 
lisions, is equivalent to a random walk with respect to 
statistical fluctuation properties. 

The stadium billiard, the Sinai billiard and the Lorentz 
gas are ergodic systems Ergodicity, i.e. long-time 
averages being equivalent to phase space averages, can 
be seen in the stadium billiard in the following way: the 
action of a trajectory for sufficiently long time (many 
collisions with the wall) can be decomposed into pieces 
corresponding to shorter trajectories between pairs of 
boundary points. This represents an average over phase 
space. Each action element follows a Gaussian distribu- 
tion and due to the strong law of large numbers, the total 
action converges against the mean value of the Gaus- 
sian. 

This work may have the following implications: (i) 
Present a new efficient method to measure global 
chaotic behavior in classical systems, (ii) Shed light on 
ergodicity /ergodicity breaking in Hamiltonian flow, (iii) 
Allow a unified description of both, quantum and classi- 
cal chaos, (iv) Help understanding why quantum chaos 
is typically weaker than classical chaos (e.g. via an ef- 
fective quantum action lll2Lll3[l ). (v) Help in searching 
for a proof of Bohigas-Giannoni-Schmit conjecture. 
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